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HIGH-ENERGY  IONS  IN  LASER-PLASMA  INTERACTIONS 

Historically,  the  high  energy  ions  in  a laser  produced  plasma  are 
defined  as  a small  group  of  ions  transporting  a significant  fraction  of 
the  absorbed  laser  energy.  Most  expansion  models1’2  suggest  that  the 
high  energy  ions  are  the  direct  consequence  of  the  expansion  of  the 
high  energy  electrons.  However,  these  models  have  only  considered  the 
acceleration  of  a single  ion  species.  Some  recent  measurements3’4  of 
the  asymptotic  energy  distribution  of  high  energy  ions  from  a multi-ion 
species  plasma  (CH2)  suggest  that  the  relative  behavior  of  the  different 
ions  species  might  lead  to  further  information  about  the  expansion  and 
the  acceleration  mechanism  of  the  ions.  In  our  model,  we  therefore 
consider  the  ambipolar  expansion  of  hot  electrons  with  two  relatively 
cold  ion  species  from  a pressure  gradient.  We  will  then  show  by 
comparison  with  experiments  that  such  a pressure  gradient  model  can 
account  for  the  energy  and  the  relative  behavior  of  these  high  energy 
ion  species. 

A typical  energy  distribution  of  the  high  energy  ions  measured 
from  a Nd-laser  target  irradiation  is  shown  in  Fig.  1.  Mutlipeak 


structure  on  the  H (hydrogen)  ion  energy  distribution  has  been  observed  «e  Section 


between  1C  and  100  keV.  ’ The  group  of  H ions  with  E/Z  (energy 
divided  by  the  charge)  higher  than  the  C (carbon)  ions  transports  as 
much  energy  as  the  whole  C ion  group.  The  presence  of  C J and  C+4 
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ions  at  almost  the  same  energy  as  the  C 0 ions  is  consistent  with  some 

atomic  processes,  such  as  recombination  or  charge  exchange,  occurring 

during  the  expansion.  We  will  therefore  assume,  in  the  model  discussed 

below,  that  all  the  C ions  were  fully  stripped  during  the  acceleration 

phase  of  the  expansion.  A detailed  discussion  of  these  observations  can 

be  found  elsewhere.3  In  this  paper,  we  will  concentrate  mostly  on  the 

•f* 

energy  partition  between  H and  C ions  and  show  that  it  is  consistent 
with  a pressure  gradient  accelerating  two  different  ion  species. 

The  1-D  plasma  is  model  ed  using  one  hot  electron  background  and 
two  relatively  cold  ion  fluids.  The  ion  density  profiles  for  our 
initial  value  problem  are  shown  in  Fig.  2a  (dashed  lines).  Both 
density  profiles  have  initially  the  same  exponential  scale  length  and 
a velocity  negligible  with  respect  to  the  final  velocities. 

The  three  species,  one  electron  and  two  ion,  are  described  by  a 
standard  set  of  conservation  equations.  Each  ion  fluid  satisfies  a 
continuity  equation 
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where  i = 1,2.  The  momentum  equation  for  the  hot  electrons  is 
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where  0 is  the  ambipolar  potential  and  p^  = ngkTe,  the  electron 
pressure.  The  momentum  equation  for  both  ion  fluids  is 
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We  also  assume  that  the  density  gradient  scale  length  l is  much  greater 
than  the  electron  Debye  length  and,  therefore,  replace  the  Poisson 
equation  for  the  ambipolar  potential  by  a quasi-neutrality  condition 

n = Zn  + Zn. 

e ii  22 

No  collisional  effects  have  been  included  in  this  model.  The 
electron-ion  collisions  can  be  neglected  due  to  the  high  electron 

temperature.  An  initial  ion  temperature5  of  a few  hundred  eV  also  makes 

_2 

the  viscosity  terra  (<*  T^2)  negligible  with  respect  to  the  electric  field 
in  Eq.  3.  6,7  The  ion  temperature,  although  high  enough  to  neglect 
viscosity,  remains  relatively  small  compared  to  the  electron  temperature, 
so  that  the  ion  pressure  can  also  be  neglected.  The  only  interaction 
between  the  two  ion  species  is  then  through  the  self-consistent 
electric  field  in  Eq.  3. 

An  assumption  about  the  electron  temperature  and  the  heat  flow  is 
required  to  close  the  moment  equations.  We  assume  a uniform  electron 

temperature  throughout  the  expansion  region,  i.e.,  the  heat  is  allowed 

iK 

to  flow  without  inhibition.  The  left  boundary  in  Fig.  2a  is  an 
impenetrable  wall  with  u^  = 0,  bpe/^x  = 0 and  no  heat  is  allowed  to  flow 
across  that  boundary.  The  total  energy  is  then  aonserved  by  reducing 
the  electron  temperature  according  to  the  rate  of  change  of  ion  kinetic 
energy. 

Equations  1 to  3 have  been  self-consistently  solved  numerically 
with  the  quasineutrality  and  energy  conservation  conditions  using  an 
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FCT  algorithm8  on  a sliding-zone  grid.  Figure  2a  shows  the 
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evolution  of  the  ion  density  profiles  for  a CH^  plasma  after  3.6  t 
(where  r is  the  density  gradient  scale  length  divided  by  the  ion  sound 
speed).  As  discussed  above,  all  the  C ions  are  assumed  to  be  completely 
ionized.  Because  of  the  more  rapid  H ion  expansion  relative  to  C 
ions  a bump  has  developed  on  the  H+  density  profile.  From  Fig.  2b-c, 
one  sees  that  this  H+  group  transports  an  important  fraction  of  the  ion 
expansion  energy  at  an  energy  about  three  times  the  initial  electron 
temperature.  In  fact,  most  of  the  ion  expansion  energy  is  contained  in 
a small  fraction  of  the  ions  with  energies  higher  than  the  initial 
electron  temperature. 

About  75^  (Fig.  2c)  of  the  electron  thermal  energy  remains  after 
3.6  t but,  because  of  the  much  weaker  density  gradient,  the  ion 
acceleration  by  the  ambipolar  electric  field  is  greatly  reduced.  The 
remaining  thermal  energy  will  then  be  dissipated  via  other  channels, 
such  as  heat  transport  toward  the  overdense  region  and  ion  heating  with 
eventual  acceleration  of  slower  ions. 

The  energy  partition  between  the  H+  and  the  C+s  ions  can  also  be 
predicted  from  the  momentum  equation  (Eq.  3).  For  both  ion  species,  the 
rate  of  change  of  ion  kinetic  energy  (E^  = £ n^v^)  is  initially  given  by 
(from  Eq.  3) 
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where  and  are  respectively  the  ion  atomic  number  and  charge  state 
and  6 is  the  ambipolar  electric  field.  Since  most  of  the  ion 
acceleration  occurs  before  the  separation  of  the  species,  the  electric 
field  and  the  ion  acceleration  time  are  essentially  the  same  for  both 
ion  species.  Therefore,  Eq.  4 yields  the  following  relation  for  the 
energy  of  the  ion  species. 

E A 

7^  (r^)  = const.  (5) 

Zi  Zi 

A different  scaling,  E^/Z^  = const. , would  be  valid  only  if  the 

different  ion  species  are  accelerated  through  the  same  electrostat ic 

potential  for  different  acceleration  times.  From  Eq.  5 a C ion  is 

+ 

then  expected  to  have  only  three  times  more  energy  than  a H ion  and 
a CH^  target  should  yield  1.5  times  more  energy  into  the  C ions  than 
the  H+  ions.  Figure  2c  shows,  however,  a C+s  to  H+  ion  energy  ratio 
close  to  1.2.  The  extra  energy  (2C^)  found  in  the  H+  ions  from  the 
simulation  is  caused  by  different  expansion  velocities  of  the  two  ion 
species.  The  H+  ions,  because  of  their  higher  velocity,  can  take 
advantage  of  the  pressure  gradient  set  up  by  the  slower  C ions. 
Nonetheless,  the  rate  equation  (Eq.  h-) , « Z2  /A,  remains  what  basically 
determines  the  energy  of  the  fast  ions. 

As  mentioned  above,  we  assume  at  3.6  t that  the  acceleration 
phase  of  the  ion  expansion  is  terminated  and  then  calculate  the 
asymptotic  energy  distribution  for  both  ion  species.  The  result  is 
shown  in  Fig.  3a  assuming  an  initial  electron  temperature  of  22.5  keV. 
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Figure  Jb  shows  the  experimental  energy  distribution  deduced  from 
Fig.  1.  Again  the  partial  current  of  C*5  and  C 4 were  assumed  to  be 
C during  the  acceleration  phase.  The  22.5  keV  initial  electron 
temperature  was  chosen  to  fit  the  H peak  from  the  expansion  model  to 
the  high  energy  H+  peak  from  the  experiment.  As  can  be  seen  in  Fig.  5b, 
the  experiment  produces  a higher  ratio  of  H+  to  C ions  than  expected 
from  a CH^  target  over  the  energy  range  considered.9  This  is,  however, 
in  qualitative  agreement  with  the  expansion  model.  The  expansion 
energy  was  shared  almost  equally  in  the  expansion  model  between  H+  and 
C ions.  Within  experimental  uncertainties  this  result  is  also 
consistent  with  the  experiment. 

A typical  ion  acceleration  time  scale  for  an  initial  electron 
temperature  of  ~ 20  keV  and  scale  length  of  a few  microns  is  a few 
tens  of  picoseconds.  The  relatively  high  electron  temperature  required 
to  explain  the  data  has  then  to  be  maintained  for  only  a short  time 
and  is  therefore  a peak  temperature.  Such  a high  temperature  is 
indeed  consistent  with  simulations  using  an  artificially  inhibited 
thermal  conductivity1  or  including  self-generated  magnetic  fields. 6,10 
A different  approach  would  have  been  to  use  a two  electron  temperature 
model  with  hot  electrons  controlling  the  ambipolar  expansion  of  the 
ions. 11 

Our  expansion  model  does  not  address  the  observation  of  the  multi- 
peak structure  of  the  H+  ion  distribution.  This  feature  could  be  due 
to  spatial  or  temporal  variation  of  the  pressure  gradient.  The  high 
energy  tail  of  the  asymptotic  high  energy  ion  distribution  can  also 


be  treated  with  an  approach  different  from  that  of  the  expansion 
model  suggested  above  (vacuum  interface  problem).12’13 

In  summary,  the  electron  pressure  gradient  model  can  account  for 
the  energy  and  the  relative  behavior  of  high  energy  ion  species.  In 
the  context  of  a multi-species  plasma  expansion  we  have  also  shown  that 
the  rate  of  change  of  ion  expansion  energy  («  Z2/A)  determines  the 
energy  partition  among  the  fast  ions.  A further  test  of  this  model 
will  then  be  to  compare  the  asymptotic  energy  distribution  from  a 
CD2  target  to  that  from  a CH2  target.  According  to  the  ^/A  scaling 
law,  the  D+  ions  should  transport  6 times  less  energy  per  ion  than  the 
C^13  ions. 

The  authors  wish  to  thank  S.  E.  Bodner,  H.  R.  Griem,  and  F,  S. 
Felber  for  their  interesting  comments  on  this  work  and  D.  G.  Colombant 
for  his  valuable  assistance  in  the  numerical  calculations. 
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E/Z  (keV) 


Fig.  1 - High-energy  ion  current  distribution  vs  energy  per  unit 
charge.  The  numbers  represent  the  Z of  C+z  ions  and  + the  H+ 
ions.  Laser  parameters:  5 x 10^  W/cm^  at  1.06  ^m  on  a planar 
CH2  target. 


Fig.  2 - Ion  densities,  ion  energies  and  energy  densities  versus 
distance  at  t = 3.  6 . The  dashed  lines  represent  the  initial  den- 

sity profiles.  The  ion  densities  are  normalized  to  the  initial  H+ 
plateau  density.  The  non-dimensional  units  are:  t - 11/ cg,  cg  = 
(kTeo/mp)  1 /2>  where  Teo  is  the  initial  electron  temperature  and 
t the  initial  density  gradient  scale  length.  The  percentages  in 
parenthesis  gives  the  energy  partition. 


0.0 


Fig.  3 - Asymptotic  high  energy  ion  distribution  (a)  from  the 
expansion  model  and  (b)  from  the  experiment.  The  energy 
scale  uses  the  A/Z^  normalizing  factor  from  the  rate  equation. 
The  error  bars  are  from  uncertainties  in  the  data  reduction 
from  current  to  dN/dE.  the  number  next  to  H+  refers  to  the 
energy  partition  among  the  high  energy  ions. 
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